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To examine thermoelectric (TE) properties of a semiconductor nanowire (NW) network, we pro¬ 
pose a theoretical approach mapping the TE network on a two-port network. In contrast to a 
conventional single-port (i.e., resistor) network model, our model allows for large scale calculations 
showing convergence of TE figure of merit, ZT , with an increasing number of junctions. Using 
this model, numerical simulations are performed for the BUTe 3 branched nanowire (BNW) and 
Cayley tree NW (CTNW) network. We find that the phonon scattering at the network junctions 
plays a dominant role in enhancing the network ZT. Specifically, disordered BNW and CTNW 
demonstrate an order of magnitude higher ZT enhancement compared to their ordered counter¬ 
parts. Formation of preferential TE pathways in CTNW makes the network effectively behave as 
its BNW counterpart. We provide formalism for simulating large scale nanowire networks hinged 
upon experimentally measurable TE parameters of a single T-junction. 


PACS numbers: 72.20.Pa, 62.23.Hj, 73.63.-b, 63.22.-m 

I. INTRODUCTION 

Environmental concerns and demands to explore alter¬ 
native energy sources have been currently motivating the 
search for new efficient thermoelectric (TE) materials 
TE efficiency is characterized by the dimensionless figure 
of merit, 

ZT = , (1) 

l^eh T Kph 

which depends on the temperature, T, electrical con¬ 
ductivity associated with electrons and holes transport 
a e h, Seebeck coefficient, S', and thermal conductivity due 
to the electrical charge carriers, n e h, and the phonons, 
Kph^- Accordingly, enhanced performance of TE mate¬ 
rials requires a high power factor, cr^S 2 , and at the 
same time suppressed thermal conductivityP Nanostruc¬ 
turing is a promising way to achieve this goal. Reduc¬ 
tion in dimensionality of the nanomaterials can be ex¬ 
ploited to quantize their electronic structure (e.g., re¬ 
sulting in density of states singularities) in order to en¬ 
hance the power factor. 4 The surface and interfaces fa¬ 
cilitate phonon scattering, subsequently decreasing asso¬ 
ciated thermal conductivity. 4 ® 

Semiconductor nanowires (NWs) are suitable can¬ 
didates for reaping the TE benefits of reduced 
dimensionality.^ Reduced dimensionality of NWs dra¬ 
matically changes their electronic density of states which 
in turn was expected to improve the power factor a e hS 2 . 

In fact, subsequent studies demonstrated that improve¬ 
ment in ZT can rather be achieved by suppressing 
the phonon contribution to thermal conductivity due to 
boundary scattering effects, surface roughness, and lat¬ 
tice anisotropyP Specially developed NW synthesis and 
characterization techniques along with theoretical model¬ 
ing have been applied to a broad range of semiconductor 


materials from groups IV, IY-VI, V-VI, II-VI, III-V to 
examine their TE capabilities. 9 13 Bulk doped Bi^-iTe^ 
is currently among the most efficient and commercially 
used TE materials with ZT ~ 1.0 at T ~ 300 KP How¬ 
ever, experimental measurements on large , ^ 50 nm, di¬ 
ameter Bi x _i Te x N Ws report a low value of ZT ~ 0.1 
at T ^ 400 KJ 14 " 15 I Theory predicts that a reduction in 
diameter to 0.5 nm should enhance Bi 2 Te 3 NW perfor¬ 
mance to ZT = 14 at T = 300 K which rapidly decreases 
to ZT ^ 1 with the diameter increase to ^ 4 nmP 
The capabilit y of gr owing thin Bi^-iTe^ NW has been 
demonstrated. 9 17 18 

Use of NWs in TE devices calls for scaling-up their 
remarkable TE properties to the macroscopic level via 
assembling them into functional blocks. Various archi¬ 
tectures for TE devices utilizing parallel arrays of NWs 
imbedded into different matrices and array growing tech¬ 
niques have been proposedP^Use of NW networks with 
a large amount of junctions instead of parallel NW ar¬ 
rays can lead to an enhancement of device performance, 
e.g., due to additional phonon scattering at the network 
junctions. For example, growth of ultra-thin, ^ 8 nm 
in diameter, Bi 2 Te 3 NWs with their subsequent fabri¬ 
cation into bulk pellets using a spark plasma sintering 
method has been reported.^ The TE characterization of 
these pellets shows an enhanced peak ZT ~ 0.96 at room 
temperature. Recently a group of researchers at Sandia 
reported creation of TE nanowire arrays with uniform 
composition along the length of Bi^-iSb^ nanowires and 
across the spread of the NW array, which potentially can 
include hundreds of millions of NWsP 

Optimization of TE devices based on NW networks re¬ 
quires modeling approaches capable to account for scal¬ 
ing up of TE properties of individual NWs to the macro¬ 
scopic device level. The initial step in bridging these 
well separated length-scales is to consider the TE prop- 
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erties of the mesoscopic NW assemblies, which can sub¬ 
sequently be used as the building blocks of macroscopic 
devices. The large number of junctions in branched NW 
mesostructures makes full scale atomistic calculations not 
feasible, thus calling for approximate theoretical mod¬ 
els. In one example, atomistic calculations were com¬ 
bined with an averaging procedure to examine the ef¬ 
fect of a single ultra-thin (< 2 nm in diameter) Si NW 
decoration with either Si branches or chemically bound 
alkyl molecules.^ 1 This study demonstrated that for 100 
branches, ZT acquires an almost four fold increase as 
a result of enhanced phonon scattering. To address the 
scaling issues in modeling the NW networks in Sec. [TI| 
we propose a generalized approach based on mapping a 
NW network onto a corresponding network of two-port 
TE elements and further taking advant age of electrical 
circuitry analysis approach! : SHU] j n g ec jjj ^he method 
is applied to examine the TE properties of model net¬ 
works made of doped B^Tes NWs. Sec. IV concludes 
our discussion. 


II. NW NETWORK MODEL 


Let us consider a TE NW network. For the tempera¬ 
ture range of 250 — 400 K used in the calculations below 
both electrical charge carriers and thermal transport oc¬ 
cur in the diffusive regime leading to small local varia¬ 
tions of electrochemical potential and temperature across 
the network. This allows us to partition the network into 
segments each of a size larger than the charge carriers and 
phonon mean free paths. 36 Each segment is terminated 
by hot and cold junctions with other segments and/or 
electrical contact as shown in Fig. [l] (a). We assign the 
difference in electrochemical potential, AV = Vh — Vc, 
and temperature, AT = Th — Tq to each segment result¬ 
ing in electric, J e h, and thermal, Jq. currents. In this 
case, the linear response theory yields^! 

J = YV (2) 


as the TE two-port impedance! 37 ! In such a network any 
two segments connected in series satisfy the sum rule for 
the impedance matrix 

Z = Zi + Z 2 (6) 

and the segments connected in parallel satisfy the sum 
rule for the admittance matrix 

Y = Yi + Y 2 . (7) 

The TE admittance matrix can be represented in terms 
of transport coefficients, for electrons, Tf and holes, T/, 
where i = 0,1, 2P^ By defining Li = Lf + T/, one can 
write the admittance matrix as 

-y _ (e 2 feh{ri)L o ef eh (r[)L\ \ . > 

\efeh(v) L i feh(v)L 2 + f P h(t,ri)KphTj' 

where, e denotes the absolute value of the electron charge, 
T = (Th+Tc )/2 is the averaged temperature, and f e h(v ) 
(fphiViO) electron-hole (phonon) reduction func¬ 

tion whose physical meaning is discussed below. 

For a single NW, the electron/hole transport coeffi¬ 
cients entering Eq. (|8| can be represented in the following 
general form 


L’(n,T) = l J T j {E){±E T lA i [jPjdE, (9) 

— OO 

where i = 0,1,2 and j = e, h. Tj(E) is the energy de¬ 
pendent transmission function and f(E, (i) = {exp[(T — 
/Lj/ksT] + l} -1 is the Fermi-Dirac distribution depend¬ 
ing on the averaged chemical potential, /q and averaged 
thermal energy, Here, the upper sign stands for 

the conduction band electrons (j = e), and the lower 
one for the valance band holes (j = h). The prefac¬ 
tor 2 accounts for the spin degeneracy. The transport 
coefficients are connected to the electrical conductivity, 
electron/hole contributio n to thermal conductivity, and 
the Seebeck coefficient as 27 28 


where the column vector of the TE current is 



and the TE potential is 


V = 


\AT/TJ ' 


( 3 ) 


( 4 ) 


Taking into account the two-component nature of these 
vectors, we propose to use the circuit analysis approach 
by treating segments of a network as TE two-port ele¬ 
ments. This allows us to identify the matrix Y as the 
TE two-port admittance, and its inverse 

Z = Y -1 (5) 


2 T 

Oeh — ^ Lo, 

L‘2 l\ 


K eh rjn 


TLn 


S = 


L\ 


eTL o’ 


( 10 ) 

( 11 ) 

( 12 ) 


respectively. Finally, the phonon conductivity in Eq. 
can be expressed as 


Kph — 


h 2 


2irkBT 


OO 

J d0JL0 2 T ph {uj) 


exp (huj/ksT) 
[exp ( huj/ksT ) — 1] 


2 ’ 


(13) 


where T p h{w) is the phonon transmission function. 

Next, we discuss the reduction functions fehiv) an d 
fph{v , 0 entering Eq. (|8| that account for the junction 
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FIG. 1: (a) NW segment terminated by hot and cold junctions, (b) The segment representation as three serial connected two- 

port elements, (c) BNW of N — 1 RUs and a single NW segment at site n — 1. Dotted boxes mark BNW portions described by 
the admittance matrices Yi,..., Y n _i, Y n . (d) Associated representation of n-th RU in terms of two segments from panel (a) 
with the admittances Y n ,i and Y n; 2 - (e) CTNW of N generations each containing 2 9 ~ x RUs, where g is the generation index. 

Dotted boxes mark the network portions each characterized by a set of admittance matrices, i.e., {YiV,4n=i , {Y 9 +l,4n=l, 
and {Y g , n }n 9 =i . (f) n-th RU from generation g modeled by two NW segments from panel (a) with the admittance matrices 
Yg jnj i and Y^ jnj 2 * 


induced drop in the electrical and thermal conductivities, 
respectively. Reduction in the electron/hole and phonon 
transmission functions associated with a single NW junc¬ 
tion results in a drop of the electrical and thermal con¬ 
ductivities. To account for this effect we introduce as¬ 
sociated single junction transmission functions using the 
following equalities 

tj(E) = (1-47 'j(E) (14) 

% h {u) = (1 -t)T ph {u), (15) 


where 0 < g < 1 and 0 < £ < 1 are the electron/hole 
and phonon reduction parameters, respectively. Tj (' T p h ) 
is electron/hole (phonon) transmission function of NWs 
forming the junction. Substitution of the junction trans¬ 
mission functions into Eqs. (|9| and (13) and comparison 
with Eq. (J8| allows one to identify single junction reduc¬ 
tion functions as f e h(v) = 1 — V and f P h(v) = 1 — £• 
Details on the reduction parameters determination are 


provided in Sec. Ill Note that the junction reduction 
parameters and associated reduction functions are as¬ 
sumed to be energy independent. This approximation 
results in the junction independent values of the Seebeck 
coefficient.^ 


Use of serial and parallel connection rules (Eqs. © and 
([7])) for a NW network results in a complex dependance of 
the reduction functions entering the admittance matrix 
(Eq. (|8|). In particular, the phonon reduction function 
becomes dependent on both electron/hole and phonon 
reduction parameters as indicated in Eq. (|8|. To model 
NW networks, we introduce a network segment (denoted 
by a set of indices s) represented by a homogenous NW 
truncated by two junctions as shown in Fig. [l] (a). Such 
a segment maps on a serial connected three two-port ele¬ 
ments shown in Fig.[l](b). The central element describes 
the homogeneous part of the NW with f e h = f p h = 1 
while the right and left elements account for an identical 
admittance drop at the junctions. The associated reduc- 
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tion functions are f eh {Vs) = 1 ~Vs and f ph (Vs^s) = 1 -&■ 
Implementing the serial connection rule (Eq. ([6])), we find 
the following expressions for the segment reduction func¬ 
tions 


fehiVs) 


fph{£.sj Vs) 


1 _ 2 4 rieML* - LqL*) 

3 9^ s 9 Lj — Lq(L2 + e 2 K p hT)' 1 


The obtained expression is exact. Starting with the ini¬ 
tial condition Yi = Yqi and performing n = 2,... N 
subsequent iterations, one can calculate the total BNW 
admittance Y bnw = Y^. Subsequently, the BNW 
figure of merit can be estimated by substituting Y = 
Y bnw into Eq. (18). 


B. Cayley tree NW network 


Note that retaining only the linear term in Eq. yields 
an identical functional form for f e h and f p h . This makes 
the electron/hole and phonon TE channels uncoupled 
and the network can be modeled using a conventional 
single-port element (i.e., resistor network) meth od. 11 
However, in the numerical calculations (Sec. |III| ), the 
cross-term is retained to account for the mixed elec¬ 
tron/hole and phonon TE channels. 

Using Eqs. (10)-(|l2|) one can easily find that the TE 
figure of merit (Eq. TO) can be expressed in terms of the 
TE admittance matrix as 


ZT = 


V 2 

*12 

detY 


(18) 


where Y\2 denotes the off-diagonal matrix element of the 
admittance matrix given by Eq. (|8|. Calculation of the 
figure of merit requires knowledge of the total network 
admittance matrix which is dependent on the network 
structure. Next we derive expressions allowing us to eval¬ 
uate the admittance matrix for two model systems such 
as the branched NW (BNW) and binary Cayley tree NW 
(CTNW) network. 


A. Branched NW 


The BNW, shown in Fig. [l] (c), contains T-shaped 
repeat units (RU) denoted by the index n = 2, ...,7V 
and a single NW segment at the site n = 1. All the 
cold contacts (blue dots) are grounded. The admittance 
of BNW portion containing n — 1 RUs (dotted box) is 
denoted by Y n _i. According to panel (d), each RU 
contains two segments denoted by the index m = 1,2. 
Each segment is characterized by the admittance matrix 
Y nm = Y (rjnm,€nm) with the reduction functions given 
in Eqs. (16) and where we set s = {n, m}. By taking 
into account that segment m = 2 of the n-th RU has a 
parallel connection with the rest of the network contain¬ 
ing n — 1 RUs and segment m = 1 is further attached via 
serial connection, we use Eqs. 0 0 to express the ad¬ 
mittance of the BNW portion containing n RUs, Y n , in 
terms of the admittance of the n — 1 RUs, Y n _i and the 
admittance of the segments forming the n-th RU, Y n?m 
(m = 1, 2) as 


Y n = 


(( 


Y n _! + Y 1 


1 , 2 ) 


r-1 

■ 71,1 


-1 


(19) 


The CTNW shown in Fig(T] ( e) contains N generations 
each denoted by the index g. Each generation is formed 
from A-shaped RUs denoted by index n = 1,... ,2^ -1 . 
For the network portion containing generations g + 1 
through N (dotted box), we introduce a set of admit¬ 
tance matrices {Y^ + i jn } with n = 1...2^. According 
to panel (f), n-th RU of generation g contains two seg¬ 
ments, denoted by the index m = 1,2. Each segment is 
characterized by the admittance Y gnrn = Y ( g Qnrn , ^ qnm ), 


where the reduction functions are given by Eqs. (16) and 

jrr 


in which we identify s = {g,n,m}. By taking into 
account that both segments in each RU of generation g 
have serial connection with the attached portion of the 
network (i.e., generations g + 1 through A), we take ad¬ 
vantage of Eqs. 0-0 to express the admittance matrix 


Yg, n as 


Y — 

x 9,n — 


(y-^ + y ;l h2n _ 1 ) 
(y~1 2 + y -L +1 ) _1 


-1 


( 20 ) 


Similar to BNW, the obtained expression is exact and 
can be used for iterative calculations of the network ad¬ 
mittance. Namely, one starts with the initial condi¬ 
tions associated with the bottom of the tree YAr, n = 
(Yjv,n,i + Y7v, n ,2) where n = 1,..., 2 Ar ~ 1 and further 
proceeds to the top by applying Eq. ( |2Q| ) for generations 
g = N — 1 through g = 1. Eventually, the total CTNW 
network admittance is identified as Yctnw = Yu and 
the CTNW figure of merit is estimated by substituting 
Y = Yctnw into Eq. (18). 


III. NUMERICAL RESULTS AND DISCUSSION 

The numerical simulations are performed for the net¬ 
works made of well parameterized n-doped Bi 2 Te 3 NWs. 
Each network segment is modeled by a rectangular NW 
elongated in the [015] crystalline direction with the cross- 
section of 8 x 8 nm 2 . The electronic structure is calcu¬ 
lated using a six-valance-band and six-conduction-band 
effective mass model accounting for the band-gap tem¬ 
perature dependence. Cold and hot reservoirs are mod¬ 
eled as infinite Bi 2 Te 3 NWs. The transport coefficients 
Lf, A- 1 are evaluated using a constant relaxation time 
model. The calculations of thermal conductivity, 
accounts for the NW boundary specularity. Details of 
these calculations are given in Appendix |A| 
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FIG. 2: ZT as a function of chemical potential calculated for 
(a) the ordered CTNW of N — 5 generations and (b) the 
disordered CTNW with N = 100 generations (thick lines). 
Thin lines in both panels present corresponding behavior 
of single NW segment, (c) ZT enhancement factor calcu¬ 
lated at T=300 K using two-port (TP) network model for 
ordered/disordered BNW as a function of number of junc¬ 
tions (i.e., RUs—1) and for the ordered/disordered CTNW 
as a function of number of generations. ZT vs. number of 
junctions in the disordered BNW calculated using single-port 
network model. 


The reduction parameters have a complex dependence 
on the junctions structure (e.g., lattice mismatch in the 
branching points). Lack of necessary microscopic param¬ 
eters does not allow for their direct evaluation. Our es¬ 
timate of the transmission function for Si NW junctions 
yields the electron/hole reduction parameter r] s ~ 0.01 
and phonon thermal conductivity reduction parameter 
£ s ~ 0.15f 39 Accordingly, branching strongly suppresses 
the phonon assisted thermal conductivity and keeps the 
electrical conductivity almost intact. We adopt the 
same order of magnitude for the reduction parameters 
in Bi 2 Te 3 NW networks as the starting point in our sim¬ 
ulations. 

Below, the following two kinds of BNW and CTNW 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


& 

FIG. 3: The enhancement factor for ordered/disordered BNW 
and CTNW containing 100 junctions and generations, respec¬ 
tively, as a function of phonon reduction parameter £ s . For 
disordered BNW and CTNW the mean value of is varied 
while the standard deviation stays 10%. The electron reduc¬ 
tion parameter is fixed at r/ s = 0.05 and T = 300 K. 

are examined: (i) Ordered BNW/CTNW in which all 
junctions have identical reduction parameters g s = 0.05 
and £ s = 0.20. (ii) Disordered BNW/CTNW in which 
each junction reduction parameter is initialized according 
to the Gaussian distribution with the mean values fj s = 
0.05 and £ s = 0.20 and 10% standard deviation. 

Figure [2] (a) presents ZT dependence on the chemical 
potential, \i = (Vh+Vc)/ 2, for the ordered CTNW (thick 
lines) containing N = 5 generations. As the chemical po¬ 
tential and associated carriers concentration increase, the 
Seebeck coefficient drops while the conductivity grows re¬ 
sulting in a well known ZT maximum. 2 A similar plot for 
the disordered CTNW of N = 100 generations is shown 
in panel (b). ZT variation for single NW segment is plot¬ 
ted in panels (a) and (b) by thin lines for comparison. 
Introducing the peak enhancement factor as 

max ( ZT bnw /ctnw ) — max (ZT^w) 

max (ZTnw) 
u 

one finds that that branching produces about 2% and 
20% enhancement in ZT peak values for the ordered and 
disordered CTNWs, respectively. Ordered and disor¬ 
dered BNWs with 5 and 100 junctions have absolutely 
the same ZT behavior as their CTNW counterparts. 

In Fig. [ 2 ] (c), we plot ZT enhancement factor for 
the ordered/disordered BNW and CTNW depending on 
the number of junctions and generations , respectively. 
The ordered BNW/CTNW demonstrate identical ZT en¬ 
hancement for all values of junctions/generations and 
saturate to ^ 2% enhancement value merely at 5 junc¬ 
tions/generations. The disordered BNW/CTNW reach 
close ZT values at around five junctions/generation and 
continue monotonous growth together. Our calculations 
(not shown in the plot) provide saturation of the en¬ 
hancement factor to ^ 20% after the number of junc¬ 
tions/generations exceeds 100. 
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FIG. 4: Temperature vs chemical potential phase di¬ 
gram. Lines separating areas marked with different shades 
of blue represent constant values of the parameter a = 
m.in{eLi/ e 2 Lp , eLi/ k p hT, L 2 / k p hT} reflecting how well in¬ 
equalities (|Bl|)-( B3) are satisfied. 


Identical variation of ZT values for BNW and CTNW 
with an increasing number of junctions and generations, 
respectively, suggests that within CTNW there is a pref¬ 
erential pathway for electrical and thermal transport. 
This pathway can effectively be considered as BNW 
whose number of junctions naturally coincides with the 
number of generations. By taking into account that re¬ 
alistic networks of NW have a more complex structure 
than adopted model network, one can develop a strategy 
of modeling realistic networks via finding the preferential 
pathways for TE current and mapping them on BNWs. 

Use of the NW networks in TE devices requires sig¬ 
nificant enhancement of ZT. Therefore, in Fig. [3] we 
show how the enhancement factor for 100-generation 
(100-junction) ordered and disordered CTNW (BNW) 
changes with the increase of the phonon reduction pa¬ 
rameter. While for the disordered network this value can 
reach 150% and higher, the ordered network shows no 
more than 20% enhancement. The disorder effect is a 
reduction of effective pathways for charge carriers and 
phonon transport. Since, the phonon reduction parame¬ 
ter exceeds the carriers one by an order of magnitude, the 
nonlinear dependance of the phonon reduction function 
on this parameter amplifies the thermal current reduc¬ 
tion and subsequently results in the drop of thermal con¬ 
ductivity. This clearly demonstrates that the disordered 
network outperforms its ordered counterpart by about an 
order of magnitude in ZT enhancement and provides a 
good candidate for device applications. 

Transport properties of TE networks and devices are 
typically modeled by considering the electrical and ther¬ 
mal currents independently. This independent channel 
regime allows for mapping a TE network on electrical and 
thermal single-port (i.e., resistor) independent networks. 
In particular, this approach has been used to examine ZT 
growths with an increased number of junctions for deco¬ 
rated Si NWs. 21 Boundaries of the independent channel 
limit for our model are formulated in Appendix [B] and 


FIG. 5: Temperature vs chemical potential phase digram. 
Dashed lines show constant values of the ratio n e h/TS 2 a e h- 


plotted in Fig. [4] for adopted Bi 2 Te 3 NW parameters. 
Blue shaded area on the phase diagram covers regions 
where the approximation holds. The smaller value of the 
scaling parameter <a, the better approximation works. 

To compare our two-port model with a conventional 
single-port one, Fig. [ 2 ] (c) presents ZT enhancement fac¬ 
tor vs. the number of junctions in the disordered Bi 2 Te 3 
BNW network calculated using both models. The plot 
shows that both models yield the same ZT behavior for a 
small (about 5) number of junctions. After that point the 
single-port network model shows much steeper growth 
than the two-port one. Our analysis reveals that in the 
limit of an infinite number of junctions, ZT calculated 
with the help of single-port network diverges. Thus we 
argue that the single-port network model is a good ap¬ 
proximation for small networks whereas modeling of ex¬ 
tended mesoscale networks should be addressed using the 
two-port network approach. 

Finally, accurate modeling of TE NW networks re¬ 
quires precise knowledge of the junction reduction pa¬ 
rameters. In the uncoupled TE channel regime (resistor 
networks) r] s and can be directly related to the single 
junction impedances, Z e h and Z p h (see Appendix |b| for 
details) as 


V = 

1 Z e h 

Zeh ’ 

(21) 

Z = 

1 _ Z p h 

Zph ’ 

(22) 


where Z e h and Z p h denote homogenous NW impedances. 
In principle, these quantities can be determined experi¬ 
mentally by measuring transport coefficients entering the 
impedance matrices in the uncoupled channels regime. 
However, according to Fig. [4j values of the chemical po¬ 
tential of (i ~ 200 /ieV and the temperature ranging be¬ 
tween T=250-350K corresponding to the maximum ZT 
do not correspond to the desired uncoupled regime. On 
the other hand, the value of Eqs. (21) and (22) is that 
they provide a simple physical interpretation of the in¬ 
troduced reduction parameters. 
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In a general case of two-port network, the reduc¬ 
tion parameters can be connected to experimentally ob¬ 
tained values of electrical and thermal conductivities for 
equal length single NW and a three segment T-junction 
(Fig. 0 (d), N = 2). As derived in Appendix [Cj the 
junction reduction parameters assume the simple form 


T) 

i 


^{^eh Oeh) 
eh &eh 
hjph) 


Kph 


(23) 

(24) 


where K p h (k p h) and a e h (cr e h) are measured single NW 
(the NW T-junction) thermal and electrical conductiv¬ 
ities, respectively. Use of the prop osed expressions as¬ 
sumes that K e h < TS 2 cr e h (Eq. ( |C5| )). Such a regime can 
be reached for a specific range of the chemical potential 
and temperature shown as the blue region in Fig. [5] 


IV. CONCLUSION 

We have developed a two-port network model to de¬ 
scribe the TE properties of the NW networks. We ar¬ 
gue that our model provides better approximation for 
mesoscale networks containing large amount of elements 
compared to conventional single-port (resistor) network 
model. Simulations based on this model are performed 
using parameters for B^Tes NWs forming BNW and 
CTNW network. They demonstrate that the disordered 
networks show about an order of magnitude higher ZT 
enhancement compared to their ordered counterparts al¬ 
though it requires extended number of generations (^ 
100 to our estimate). Specifically for considered dis¬ 
ordered BNW and CTNW, calculated ZT values show 
enhancement in the range of 20 — 150% compared to a 
single NW segment. Observed similarity in ZT growth 
with the number of junctions and generations in BNW 
and CTNW, respectively. This suggests that preferential 
transport pathways form in CTNW effectively acting as 
BNWs. This provides an approach for modeling more 
complex network structures by mapping them on a set 
of BNWs. In general, use of efficient electrical circuitry 
simulation tools combined with our two-port TE model 
should allow for direct examination of complex NW net¬ 
work structures. Simulations of junction reduction pa¬ 
rameters for complex material structures is limited. As 
an alternative, we provide an approach for their determi¬ 
nation from experimental measurements on single junc¬ 
tions in a single-port (resistor) network limit and using 
a T-junction in a general case of the two-port network 
model. 
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Appendix A: Microscopic expressions of TE 
admittance matrix elements 


In this appendix, we present microscopic expressions 
for the transport coefficients entering TE admittance ma¬ 
trix given by Eq. These coefficients describe the 

charge carrier contributions to the electrical and thermal 
conductivity and the Seebeck coefficient. Representation 
for the phonon-assisted thermal conductivity entering the 
TE admittance matrix is provided as well. Finally, we 
discuss parameterization of B^Tes NW. 

In the constant relaxation time approximation ( valid 
at T > 200 K), the Drude conductivity model yields 28 29 

li = cjY,?- 1 / 2 / 4 , ( A1 ) 

a, (3 

L 3 i = j3.Fi/2/4> “F i/2^ a g j 5 (A 2 ) 

cx.,/3 

Li = CjUbT) 2 ^ 5J- 3/2 ^t6^J-i/2^ (A3) 

a,/3 

+ U p ) 2 T- 1/2 ^ p h 


where the integer indices cq (3 stay for the quantum num¬ 
bers associated with the size quantization across the NW. 
The prefactor Cj = (! 2ejljm* y /'KAh ) yJ2m* y kBT con¬ 
tains the NW cross-section, A, carriers effective mass 
along the NW axis (identified as y- axis) , and the car¬ 

riers mobility jlj . The latter quantity is evaluated accord¬ 


ing to the Matthiesson’s rule, 1 /jlj = (l + (A/^/A/7r)(l — 
p)/( 1 + p))/jtfj, using as parameters bulk carrier mobil¬ 
ity, /A, mean free path A, and surface specularity pP^In 
Eqs. ( Al|)-(|A3|), 


/ T 171 dT 

exp(x =F W + 1 
o 

is the Fermi-Dirac fractional integral whose argument is 
the reduced chemical potential, y 3 a p = (y— T^ /3 (0))//cbT, 
depending on the carriers energy at the band edge of as¬ 
sociated dispersion curve, E° a ^(0) = E° a ^{k = 0), char¬ 
acterized by the carriers wave vector, /c, along the NW 
direction. Finally, summation over all dege nerat e elec¬ 
tron/hole pockets must be carried out in Eqs. 

To model the electronic structure of B^Tes NW, we 
adopted the effective mass k • p approach. Specifi¬ 
cally, rectangular BUTes NW grown along [015] crystal¬ 
lographic direction is considered. 18 As mentioned above, 
the NW axis is aligned along the ^/-direction so that the 
(x, 2 :)-plane contains the NW cross-section. The sides of 
the rectangle are set to be identical, a z = a x m 8 nm. 












According to this model the electron/hole energy at the 
band edge is 


E^ = ±E g (T)/2± 


2^2 


h 2 7T 


Of 


P 2 


rnf-a, 


J,x^x 


m- a 


(A4) 


'3,z z 


Here, zero energy is set at the middle of the bulk band 
gap. The temperature dependent values of the effective 
masses m* ?n with j = e, h and n = x,y,z are adopted 
from the literature! 29 -! The band gap temperature depen¬ 
dence is E g (T) = 150.0 - 0.0947T 2 /(T + 122.5) meVP 
Finally, the bulk mobilit ies are s et to fl b e = 652.8 cm 2 /Vs 
and jl h h = 345.9 cm 2 /Vs P 9131 1 321 and scale with the tem¬ 
perature as jl b e ~ T~ 17 and jd b h ~ T _2 0 P^The adopted 
specularity parameter is p = 0.95 and the mean free path 
of the charge carriers is set to A « 40 nm. 15 

The boundary roughness described by the specularity 
parameter p is crucial for the phonon transport in B^Tes 
NW. 15 To account for the specularity effect, solution of 
the Boltzmann equation is complimented with the corre¬ 
sponding boundary condition. This along with the Drude 
conductivity in place of T p h{w) results irP^ 

_ f k B \ 3 k B T 3 A ph n n( ^ ( 

Kph - \ h) 27 t 2 v 2 (1 (A5) 

0 d/T 

f e 4 exp (e)de 
J (exp(e) — l) 2 5 

where the specularity function is 



Here v is the average value of the acoustic phonon group 
velocity weighted by the phonon population factor. A p h 
is the phonon mean free path, e = kuj/k B T is normal¬ 
ized phonon energy, and ©£> is the Debye temperature. 
For the narrow wires (A p h < y/A where A is the NW 
cross-section), the group velocity and the mean free path 
become functions of e due to quantization of the phonon 
energy. However here, we consider the situation in which 
y/A A p h so that the size quantization of phonons be¬ 
comes negligible and the specularity effect dominants. 

In the numerical simulations of B^Tes NWs the follow¬ 
ing parameters are used: The transverse, longitudinal, 
and twisting branches of the acoustic phonons are charac¬ 
terized by their group velocities v = 8.47x 10 5 , 5.34 x 10 5 , 
and 3.95 x 10 5 cm/s, respectively. The phonon mean free 
path is set to A p h ~ 2.5 nm, 15 and the Debye tempera¬ 
ture is chosen so that K p h ~ 450/T for T > 2001T, and 

p = 0.95P . 


Appendix B: Single-port network limit 


We restrict our model to the regime of uncoupled elec¬ 
trical and phonon channels facilitating single-port (re¬ 
sistor) network approximation. This approximation is 
valued under the conditions 


eLi <C e 2 L 0 , (Bl) 

eLi « K ph T, (B2) 

L 2 « K ph T, (B3) 

simplifying the admittance matrix (Eq. <©) to the fol¬ 
lowing diagonal form 


(e 2 feh(v) L 0 0 \ 

\ 0 fph(£) K ph,T) 


(B4) 


Accordingly, the electrical and phonon impedances of a 
NW (i.e., feh(v) = fph(0 = 1) can be defined as Z eh = 
(Y- 1 )!! = (e 2 io) _1 and Z ph = (Y- 1 ) 22 = {n ph T)~\ 
respectively. For a single junction, the correspond¬ 
ing impedances become defined as = (Y- 1 )!! = 
(JfehWLo)- 1 and Z ph = (Y _1 ) 22 = (fphiO^hT)- 1 . 
Substitution of the reduction functions in the form 


feh(Vs) = 1 - Vs and f p h(Vs,ts) = 1 - 6 into these ex- 
pressions immediately results in Eqs. (21) and (22). 


Appendix C: Relationship between reduction 
parameters and electrical and thermal conductivities 
for NW T-junction 


In this Appendix, we examine a T-shaped junc¬ 
tion made of three identical NWs. For such a junc¬ 
tion, we denote the electrical conductivity o e h, thermal 
electron/hole-assisted conductivity k e h, thermal phonon- 
assisted conductivity k p h, and the Seebeck coefficients 
S. By using general relations (Eqs. (10)—(12)) between 
transport coefficients and the electrical conductivity, we 
represent the T-junction admittance matrix (Eq. as 


Y 2 = 


( cWe 


TSa eh /e 


\TSa eh /e T 2 S 2 a eh + Tk eh +Tk ph 


(Cl) 


Note that the effect of the junction is already included 
into the transport coefficients (Z^, i = 0,1,2) and no 
reduction functions are needed in Eq. (Cl). 


On the other hand, we introduce similar parameters 
for each NW forming the T-junction as 3 a e h = e 2 Lo, 
S = Li/(eTL 0 ), 3^ eh = L 2 /e 2 T - L 2 /(e 2 TL 0 ), and 
3 K p h- Here, the factor of 3 comes from the fact that each 
NW segment has three times larger conductivity than the 
whole NW. Further applying Eq. (19) for the T-junction 


with N = 2 along with the reduction functions in the 
from of Eqs. (16) and ( fl7| ), one obtains 
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Y 2 


/ 2a eh (r]-l) 

I e 2 (rf—2) 

1 2STa eh ( V -l) 2T 2 S 2 a eh (r)-l ) 
\ e(y-2) r)-2 


2STa eh (r 1 -l ) 

e ( 7 7 — 2) 

2T*(/^ph, ~l~ ^eh) ((£ 1) Kph (jj l)/^eh.) 

(^-2)/tp^ + (r7-2)K e h 


(C2) 


Comparing the admittance matrix elements in Eq.( Cl) 


with Eq.(C2), one arrives at the following expression for 
the reduction parameters in terms of the introduced elec¬ 
trical and thermal conductivities 


T] = 


£ = 2 


+ 


2 (&eh & eh) 

2 O'eh &eh 

2 Keh&eh 


‘ZhvphO'eh Kph&eh 

2 (fcph ^eh) 


(C3) 

(C4) 


Kph iPph T K*eh 2(a €ph l^eh)) 


Deriving this expression, we used the fact that the junc¬ 
tion does not cha nge the Se ebec k coefficient, i.e., S = S. 
As a result Eqs. (C3) and (C4) become independent of 


S. To uncouple the electron-assisted and phonon-assisted 
effects on thermal conductivity, we impose a mild con¬ 
straint 


^eh $ T(7 e h' 


(C5) 


This conduction elimi nate s k e h and K e y in the matrix 
element (Y 2 ) 22 of Eq. (Cl) and (C2), resp ectiv ely. Sub¬ 
sequently, one can set k e y = K e h = 0 in Eq. (C4) resulting 
in 


£ = 


2 (hiph Kph) 

‘Ihiph Kph 


(C6) 


Obtained Eqs. (C3) and (C6) are used in the main text 
as Eqs (23) and (24), respectively. 


* Electronic address: oroslyak@fordham.edu 
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